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Abstract 

We report new calculations of the cross sections for deeply inelastic neutrino- 
nucleon scattering at neutrino energies between 10 9 eV and 10 21 eV. We com- 
pare with results in the literature and assess the reliability of our predictions. 
For completeness, we briefly review the cross sections for neutrino interactions 
with atomic electrons, emphasizing the role of the TF-boson resonance in u e e 
interactions for neutrino energies in the neighborhood of 6.3 PeV. Adopting 
model predictions for extraterrestrial neutrino fluxes from active galactic nu- 
clei, gamma-ray bursters, and the collapse of topological defects, we estimate 
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event rates in large-volume water Cerenkov detectors and large-area ground 
arrays. 
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I. INTRODUCTION 



Neutrino observatories hold great promise for probing the deepest reaches of stars and 
galaxies Unlike charged particles, neutrinos arrive on a direct line from their source, 

undeflected by magnetic fields. Unlike photons, neutrinos interact weakly, so they can 
penetrate thick columns of matter. For example, the interaction length of a 1-TeV neutrino 
is about 2.5 million kilometers of water, or 250 kilotonnes/cm 2 , whereas high-energy photons 
are blocked by a few hundred grams/cm 2 . 

Ultrahigh-energy neutrinos can be detected by observing long-range muons produced 
in charged- current neutrino-nucleon interactions. To reduce the background from muons 
produced in the atmosphere, it is advantageous to site a neutrino telescope at a depth of 
several kilometers (water equivalent) or to observe upward-going muons. The reactions 
(vi, ue)N — > (£~,£ + ) + anything and (u£, V£)N — > (i/£, vg) + anything are the major sources 
of both the desired signal and the attenuation of the neutrino "beam" as it passes through 
the Earth en route to the detector. Through the past dozen years, improving knowledge 
of the partonic structure of the nucleon has made possible a series of increasingly refined 
predictions for the interaction cross sections []5| |T3|] . Over the same period, ideas about the 
flux of neutrinos from active galactic nuclei (AGNs) and other extraterrestrial sources have 
evolved considerably. The observation []T3| of neutrinos correlated with supernova SN1987A 
and the detection of solar neutrinos by observing the direction of recoil electrons from 



neutrino interactions [0,0 showed the promise of neutrino observatories for astrophysical 
studies. The detection of neutrinos produced by cosmic-ray interactions in the Earth's 
atmosphere |T7|-p0l has emerged as a tool for investigating neutrino oscillations |21-24 



Plans for neutrino observatories that will detect neutrinos that originate beyond Earth have 
matured to the point that it is now reasonable to contemplate instrumenting a volume of 
water or ice as large as 1 km 3 . p5|-|32|. The ground array of the proposed Pierre Auger 
Cosmic Ray Observatory ]33| would have an acceptance exceeding 1 km 3 sr of water for 
neutrino energies greater than 10 17 eV ||34|| . The Orbiting Wide-angle Light collectors project 
(owl) would place in Earth orbit a lens to study air showers initiated by > 10 20 -eV particles, 
including neutrinos 



In this paper we present new calculations of the cross sections for charged-current and 



neutral current interactions of neutrinos with nucleons, updating our results of Ref. |I2 
(GQRS96) to take account of new information about the parton distributions within the 
nucleon fl3"6| . In place of the CTEQ3-DIS parton distributions |37| that we adopted as our 
nominal set in GQRS96, we base our new calculations on the CTEQ4-DIS parton distri- 



butions [38|. The changes are modest, and only noticeable at the highest neutrino energies 
we consider: at E u = 10 21 eV, the new cross sections are about 25% smaller than those of 
GQRS96. We find that for neutrino energies up to 10 16 eV, all the standard sets of parton 
distribution functions yield very similar cross sections. At higher energies, the predictions 
rely on incompletely tested assumptions about the behavior of parton distributions at very 
small values of the momentum fraction x. The resulting uncertainty reaches a factor 2 ±l 
around 10 20 eV. 

We combine our new evaluations of the neutrino-nucleon cross sections with models for 
the flux of ultrahigh-energy (UHE) neutrinos to estimate event rates in neutrino observa- 
tories. We consider the diffuse flux of neutrinos from AGNs and the flux of neutrinos that 
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may accompany gamma-ray bursts, as well as neutrinos from cosmological sources such as 
the decay of topological defects formed in the early universe. We evaluate rates for upward- 
going muons produced in or beneath large underwater and ice detectors, and we compute 
rates for contained neutrino interactions in a km 3 volume. We have also estimated rates for 
the proposed ground array of the Pierre Auger Cosmic Ray Observatory. 

The detection of upward-going muons from AGNs looks feasible in the next generation 
of underground experiments with effective areas on the order of 0.1 km 2 . As the muon 
energy threshold increases above a few TeV, atmospheric neutrinos and muons become 
less important backgrounds. Downward and air-shower event rates look promising for km 3 
detectors, for a variety of models. 

In the next Section, we review the calculation of the neutrino-nucleon charged-current 
and neutral-current cross sections, their sensitivities to parton distribution functions, and 
the resulting neutrino-nucleon interaction lengths. We also give a brief account of neutrino- 
electron cross sections and interaction lengths. Ultrahigh-energy neutrino rates for a selec- 



tion of flux models appear in §111. Our summary and conclusions make up §IV. 



II. NEUTRINO-NUCLEON INTERACTIONS 

A. Inclusive cross sections 

We calculate the inclusive cross section for the reaction 

v^N — > n~ + anything, (2.1) 

TX -)- p 

where N = — - — is an isoscalar nucleon, in the renormalization group-improved parton 

model. The differential cross section is written in terms of the Bjorken scaling variables 
x = Q 2 /2Mu and y = vjE v as 



d 2 a 1G\ME V ( Ml, 2 



1 w 



dxdy 7i \Q 2 + M V/ 



xq(x,Q 2 ) + xq(x,Q 2 )(l-y) 2 , (2.2) 



where — Q 2 is the invariant momentum transfer between the incident neutrino and outgoing 
muon, v = E v — is the energy loss in the lab (target) frame, M and My/ are the nucleon 
and intermediate-boson masses, and Gp = 1.16632 x 10~ 5 GeV -2 is the Fermi constant. 
The quark distribution functions are 



o u v (x,Q 2 ) + d v (x,Q 2 ) u s (x,Q 2 ) + d s (x,Q 2 ) 



+s s (x,Q 2 ) + b s (x,Q 2 ) (2.3) 

q{x, Q ) = h c s {x, Q )+t s (x,Q ), 

where the subscripts v and s label valence and sea contributions, and u, d, c, s, t, b denote 
the distributions for various quark flavors in a proton. At the energies of interest for neutrino 
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astronomy, perturbative QCD corrections to the cross section formula (|2.2|) are insignificant. 
In the DIS factorization scheme appropriate to the CTEQ4-DIS parton distributions, the 
terms proportional to a s [^] in the NLO cross section contribute only a few percent, so we 
omit them. 

The ti sea is a negligible component of the nucleon over the Q 2 -range relevant to neutrino- 
nucleon scattering; accordingly we neglect it. At the energies of interest here, it is a sound 
kinematical simplification to treat charm and bottom quarks as massless. However, the 
threshold suppression of the b — > t transition must be taken into account. We adopt the 
standard "slow-rescaling" prescription [i0|, with m t = 175 GeV/c 2 . We have carried out 
numerical integrations using the adaptive Monte Carlo routine VEGAS ||41|| , and Gaussian 
techniques. 

The neutral-current cross section is of interest because it contributes to the attenuation 
of neutrinos as they pass through the Earth. Neutral-current vN interactions may also be 
significant backgrounds to the observation of the resonant formation process v e — *> W~. 
If it becomes possible to measure neutral-current reactions and characterize the neutrino 
energy, Carena, et al. |^2[ have shown that the neutral-current to charged-current ratio is 
an important discriminant of new physics. 

Within the electroweak theory, a calculation parallel to the one described above leads 
to the neutral-current cross section. The differential cross section for the reaction v^N — > 
v „ + anything is given by 



d 2 a 
dxdy 



G 2 F ME V 
2tt 



Q 2 + Mf 



xq"(x,Q 2 )+xq°(x,Q 2 )(l-y) 



(2.4) 



where Mz is the mass of the neutral intermediate boson. The quantities involving parton 
distribution functions are 



q°(x,Q 2 



2 ) u s (x,Q 2 ) + d s (x,Q 2 ) 



(Ll + L 2 d ) 



+ 



u s (x, Q 2 ) + d s (x, Q 2 



(Rl + Rl) + 



[s s (x, Q 2 ) + b s (x, Q 2 )}(L 2 d + R 2 ) + [c s (x, Q 2 ) + t s (x, Q 2 )}{L 2 U + R 2 U ) 



(2.5) 



q°(x,Q 2 ) 



u v (x,Q 2 ) + d v (x,Q 2 ) u s (x,Q 2 ) + d s (x,Q 2 ) 



(Rl + R 2 d ) 



+ 



u s (x, Q 2 ) + d s (x, Q 2 



(L 2 U + L 2 d ) + 



[s s {x, Q 2 ) + b s (x, Q 2 )](L 2 d + Rl) + [c s {x, Q 2 ) + t s (x, Q 2 )](L 2 U + R 2 U ) 
where the chiral couplings are 



L u — 1 — \xw L d — — 1 + \xw 



R u 



- x x w 



Ra 



\x w 



(2.6) 



(2.7) 



and Xw = sin 2 8w is the weak mixing parameter. For numerical calculations we have chosen 
xw = 0.226, consistent with recent measurements Again the top-quark sea is negligible. 
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We show the cross section for the charged-current reaction u^N — > /x _ + anything as a 
function of the neutrino energy E v in Figure [l] (thin solid line). At low energies the charged- 
current cross section o"cc rises linearly with E v . For energies exceeding about 10 4 GeV, the 
cross section is damped by the W^-boson propagator. We also show in Figure [l] the neutral- 
current cross section cxnc for the reaction v^N — ► + anything (dashed line), together with 
Otot, the sum of charged-current and neutral-current cross sections (thick solid line). For 
the range of neutrino energies of interest here, the charged-current results apply equally 
to the reaction v e N — > e~ + anything. The neutral-current cross section for the reaction 
v e N — > u e + anything is identical to ctnc depicted here. 

The CTEQ4-DIS parton distributions are somewhat less singular as x — > than the 
CTEQ3-DIS parton distributions we adopted as our nominal set in GQRS96. Specifically, 
the sea-quark distributions of the CTEQ4 set behave as 

xgI CTEQ4 ](x)ocx- a227 (2.8) 

near x = 0, whereas those of the CTEQ3 set behave as 

xgF TEQ3 ](x)ocx- a332 . (2.9) 

The gentler singularity of the CTEQ4 distributions implies a smaller cross section at the 
highest energies, where the predominant contributions to the cross section come from very 
small values of x. We show in Figure |2| the ratio of the charged-current cross sections calcu- 
lated using the CTEQ4-DIS and CTEQ3-DIS parton distributions. Up to E v « 10 7 GeV, 
the two evaluations agree within a few percent. At still higher energies, the CTEQ4 cross sec- 
tion falls below the CTEQ3 cross section. At the highest energy we consider, E v = 10 12 GeV, 
the ratio is 0.74. This is a small change. 

Similar calculations lead to the cross sections for vN scattering. We show in Figure [3] the 
neutral-current (dashed line), charged-current (thin solid line), and total (thick solid line) 
uN cross sections. At low energies, where the contributions of valence quarks predominate, 
the i>N cross sections are smaller than the corresponding vN cross sections, because of 
the familiar (1 — y) 2 behavior of the vq cross sections. Above E v m 10 6 GeV, the valence 
contribution is negligible and the vN and vN cross sections become equal. 

We collect in Tables | and the charged-current, neutral-current, and total cross sections 
for vN and vN interactions, respectively. For the angular distributions, characterized by 
the mean inelasticity parameter (y), we refer to the CTEQ3 values we presented in Tables 
1 and 2 of Ref. 0. 

For 10 16 eV<E„< 10 21 eV, the CTEQ4-DIS cross sections are given within 10% by 

/ E \ °- 363 
a cc {vN) = 5.53 x lO" 36 cm 2 (^^v" ) 

a NC (WV) = 2.31 x 10" 36 cm 2 (y^J 

, p N 0.363 

a tot (WV) = 7.84 x 10- 36 cm 2 [j^) 
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rp v 0.363 

a cc (vN) = 5.52 x lO" 36 cm 2 [y^sf) ( 2 ' 10 ) 



a NC (uN) = 2.29 x 10~ 3ti cm 2 (j^) 



0.363 

36 „™2 



P N 0.363 

^(PJV) = 7.80 x 10- 36 cm 2 (j^) 



B. Variant parton distributions 



The CDF Collaborations's suggestion [Q that the yield of jets in the reaction pp 



jet 1 + jet 2 + anything exceeds the rate expected in standard quantum chromodynamics has 
prompted a re-examination of the uncertainties of parton distributions at moderate and large 
values of x. In particular, the CTEQ Collaboration has produced a variant of the CTEQ4 
distributions, labelled CTEQ4-HJ, in which an enhanced gluon population at large values 
of x raises the predicted two-jet inclusive cross section. The increased gluon density at large 
x can affect the small-x sea-quark distributions at large values of Q 2 , so it is interesting to 
ask what difference the variant parton distributions would make for the UHE neutrino cross 
sections. 

We show in Figure § the UHE vN cross sections calculated with the CTEQ4-HJ parton 
distributions. They are quite similar to those calculated with the standard CTEQ4-DIS 
distributions. The ratio of CTEQ4-HJ to CTEQ4-DIS cross sections is displayed in Figure 
[5]. The difference is less than 5% up to 10 8 GeV, and is smaller than 15% at the highest 
energy we consider. It is of no consequence for neutrino observatories. 



C. Assessment 



How well is it possible to predict a{v^N — > i + anything) and o~{v^N — ► i/g + anything)? 
For E u < 10 16 eV, all the standard sets of parton distributions, by which we mean those 
fitted to a vast universe of data, yield very similar cross sections |T2l , |IB ,H],ri|, within the 
standard electroweak theory. For E v > 10 16 eV, cross sections are sensitive to the behavior 
of parton distributions at very small x, where there are no direct experimental constraints. 
At these high energies, different assumptions about x — ► behavior then lead to different 
cross sections. Judging from the most extreme variations we found in GQRS96, and from 
our present calculations, we conclude that at 10 20 eV, the uncertainty reaches a factor of 
2 ±x . 

New physics can, of course, modify the UHE cross sections. The contributions of su- 
perpartners have been evaluated by Carena, et al. [p2 |. Doncheski and Robinett [^5| have 



investigated leptoquark excitation. Bordes and collaborators [IB] have speculated that new 
interactions might dramatically increase the UHE vN cross sections, but Burdman, Halzen, 
and Gandhi have countered H7J that unitarity limits the growth of the UHE cross section. 
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D. Interaction lengths 



The neutrino beams produced in accelerator laboratories are purified by passage through 
several kilometers (water equivalent) of material, which absorbs any accompanying particles 
without attenuating the neutrino flux. The story is different at the ultrahigh energies of 
interest to neutrino astronomy. The rise of the charged- current and neutral-current cross 
sections with energy is mirrored in the decrease of the (water-equivalent) interaction length, 

^ int = 7r (TP \N ' t 2 " 11 ) 

where N A = 6.022 x 10 23 mob 1 = 6.022 x 10 23 cm 3 (water equivalent) is Avogadro's 
number. The energy dependence of the interaction lengths for neutrinos on nucleons is shown 
in Figure |6|. We show separately the interaction lengths for charged- current and neutral- 
current reactions, as well as the interaction length corresponding to the total (charged- 
current plus neutral-current) cross section. The same information is shown for antineutrinos 
on nucleons in Figure [7]. Above about 10 16 eV, the two sets of interaction lengths coincide. 
These results apply equally to v e N (or v e N) collisions as to v^N (or v^N) collisions. 

Over the energy range of interest for neutrino astronomy, the interactions of z/ e , z/^, and 
with electrons in the Earth can generally be neglected in comparison to interactions 
with nucleons. The case of v e e interactions is exceptional, because of the intermediate- 
boson resonance formed in the neighborhood of E T u es = M^//2m « 6.3 x 10 15 eV. The 
resonant reactions v e e — > W~ — > and v e e — > W~ — > hadrons may offer a detectable 
signal. At resonance, the reaction v e e — > W~ — > anything significantly attenuates a v e beam 
propagating through the Earth. The water-equivalent interaction lengths corresponding to 
the neutrino-electron cross sections computed in [|I^] are displayed in Figure |8|. These are 
evaluated as 

A ^ = a, e (^)(10/18)iV A ' (2J2) 

where (10/18)A^a is the number of electrons in a mole of water. 

We have reviewed current knowledge of the structure of the Earth in Ref. ||12|| . To good 
approximation, the Earth may be regarded as a spherically symmetric ball with a complex 
internal structure consisting of a dense inner and outer core and a lower mantle of medium 
density, covered by a transition zone, lid, crust, and oceans. A neutrino emerging from the 
nadir has traversed a column whose depth is 11 kilotonnes/cm 2 , or 1.1 x 10 10 cmwe. The 
Earth's diameter exceeds the charged-current interaction length of neutrinos with energy 
greater than 40 TeV. In the interval 2 x 10 6 GeV < £„ < 2 x 10 7 GeV, resonant i> e e 
scattering adds dramatically to the attenuation of electron antineutrinos. At resonance, 
the interaction length due to the reaction v e e — *> W~ — *> anything is 6 tonnes/cm 2 , or 
6 x 10 6 cmwe, or 60 kmwe. The resonance is effectively extinguished for neutrinos that 



traverse the Earth. In the estimates of event rates that follow in §|T|, we take account of 
the effect of attenuation on upward-going neutrinos. 
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III. ASTROPHYSICAL NEUTRINO FLUXES: UHE EVENT RATES 



Since the publication of the GQRS96 cross sections and event rates []I2"| |, several new 
models of the diffuse neutrino flux from active galactic nuclei (AGNs) have appeared |f48| , f49 . 
In addition, Waxman and Bahcall |50| have argued that it might be possible to detect 
neutrinos associated with gamma-ray bursts. New models have also been put forward for 
the production of neutrinos in the decays of generic heavy particles |51] and in the collapse 
of topological defects |)2| . 



It is therefore timely to reconsider the event rates to be expected in large-volume detec- 
tors. We focus on the production of upward-going muons in the charged-current reactions 
(v^jV^N — > + anything. Upward-going muons are free of background from the 

flux of muons produced by cosmic-ray interactions in the atmosphere. It is in any case 
advantageous to site a detector beneath several kmwe to shield it from the (downward) 
rain of atmospheric muons ||53[ . Even at 3 kmwe underground, a detector still sees more 
than 200 vertical muons km -2 s~ x sr -1 , though most of these muons are quite soft. If we 
impose the requirement that E™ in > (10 3 , 10 4 , 10 5 ) GeV, the flux is (7,3 x 10~ 2 ,6 x 10" 5 ) 



muons km - s sr [ 5i| . As the incident zenith angle of the atmospheric muons increases, 
the background flux decreases. For horizontal incidence and below, the muon rate observed 
underground should be largely background-free. There is another important reason for look- 
ing down: The few-km range of UHE muons means that large- volume detectors can observe 
charged-current events that occur not only within the instrumented volume, but also in the 
rock or water underlying the detector. Accordingly, the effective volume of a detector may 
be considerably larger than the instrumented volume, for upward-going muons. For energies 
above 40 TeV, the Earth's diameter exceeds the interaction length of neutrinos. At these 
energies it is beneficial to look for events induced by downward and horizontal neutrino 



conversions to muons 57 



A. Sources of UHE Neutrinos 

We display in Figure |9] the neutrino fluxes calculated in a number of models. Neutrinos 
produced by cosmic-ray interactions in Earth's atmosphere dominate over other neutrino 
sources at energies below a few TeV. In this energy regime, the flux of atmospheric neutrinos 
is derived from the decay of charged pions and kaons. The dotted curve in Figure shows 



the angle-averaged atmospheric (ATM) + flux calculated by Volkova p8||59l , which we 
parametrize as 

^%±^ = 7.8 x 10- 11 ( -l^y 3 ' 6 cm" 2 s- 1 sr" 1 GeV" 1 . (3.1) 
dE v VI TeV/ v ; 

The prompt neutrino flux from charm production in the atmosphere, a small additional 
component that appears above E v « 1 TeV, has been estimated recently by Pasquali, et al. 



Active galactic nuclei (AGNs) are the most powerful radiation sources in the universe, 
with luminosities on the order of 10 45±3 erg/s. They are cosmic accelerators powered by the 
gravitational energy of matter falling in upon a supermassive black hole. Protons accelerated 
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to very high energies within an AGN may interact with matter in the accretion disk, or with 
ultraviolet photons in the bright jets along the rotation axis. Charged pions produced in 
the resulting pp or pry collisions decay into muons and muon neutrinos. The subsequent 
muon decays yield additional muon neutrinos and electron neutrinos. Neutrino emission 
from AGNs may constitute the dominant diffuse flux at energies above a few TeV. 

The solid lines in Figure || show the + fluxes predicted in several contemporary 
models for neutrino production in AGNs: Protheroe's model E8[] of neutrinos produced in 



P7 interactions, denoted AGN-P96 (jry); and Mannheim's model |^9|] of neutrino production 
in pry interactions, denoted AGN-M95 (pj). The neutrino fluxes (AGN-SS91) based on 
the pioneering work of Stecker and collaborators are still consistent with measurements. We 
retain the Stecker-Salamon flux considered in GQRS96 as a baseline to indicate the changes 
in event rates due to our new evaluation of the cross section. 



A mechanism for gamma-ray bursts put forward by Waxman and Bahcall |50[ also yields 
UHE neutrinos. The isotropic flux (GRB-WB) that results from a summation over sources 
is given by 



dN ^+*» _ M ( J^-) n cm -2 B-i bt-i G eV-\ (3.2) 



dE v VI GeV 7 

with (Af = 4.0 x 10~ 13 ,n = 1) for E v < 10 5 GeV, and (X = 4.0 x 10" 8 ,n = 2) for 
E u > 10 5 GeV. Sigl, et al. [51j have explored a class of models of exotic heavy particle 



decays that ultimately lead to neutrinos in the final state. These are referred to as top-down 
(TD-SLSC) models. In particular, we consider the model in which the heavy X-particles 
have mass Mx = 2 x 10 16 GeV/c 2 , which may arise from the collapse of networks of ordinary 
cosmic strings or from annihilations of magnetic monopoles. An interesting feature of this 
model is that the highest energy cosmic rays are photons. In the conventional topological 
defects model, the network of long strings loses its energy to the gravitation radiation. 



Wichoski, et al. |52[ have proposed a model in which particle production is the dominant 
channel through which energy is lost. Even with the observational limits from Frejus and 
Fly's Eye as constraints, this non-scaling model (denoted TD-WMB) produces a much higher 
neutrino flux than the TD-SLSC model of Sigl, et al, for the largest possible string mass 
density. The fluxes predicted by these exotic sources are shown as dashed curves in Figure 

The Super-Kamiokande Collaboration has recently presented evidence that muon neu- 
trinos produced in the atmosphere oscillate into tau neutrinos or sterile neutrinos during 
their passage through the Earth |JTJ . At the neutrino energies of interest for the detection of 



extraterrestrial sources, E v > 1 TeV, the probability for atmospheric neutrinos to oscillate 
en route to the detector, 

■> ( fi A™ 2 L 1 TeV\ , s 

is less than 10 -3 for a path length L comparable to an Earth diameter (w 13, 000 km), if the 
neutrino mass-squared difference Am 2 rs 2.2 x 10~ 3 eV 2 . For UHE neutrinos from distant 
sources, the oscillation probability may become interestingly large. 
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B. Neutrino nucleon interactions 



The upward-muon event rate depends on the v^N cross section in two ways: through 
the interaction length that governs the attenuation of the neutrino flux due to interactions 
in the Earth, and through the probability that the neutrino converts to a muon energetic 
enough to arrive at the detector with larger than the threshold energy E™ in . 

The probability that a muon produced in a charged-current interaction arrives in a 
detector with an energy above the muon energy threshold E™ m is given by 

PriE,,, Ef») = N A a cc {E v )(R{E v] Ef»)), (3.4) 



where (R(E U ; E™ m )} is the average range of a muon in rock [55] and N A is Avogadro's 
number. Although the Earth is transparent to low-energy neutrinos, an Earth diameter 
(1.1 x 10 10 cmwe) exceeds the interaction length of neutrinos with energies higher than 



about 40 TeV. For the isotropic fluxes presented in § [III A| , we represent the attenuation of 



neutrinos traversing the Earth by a shadow factor that is equivalent to the effective solid 
angle for upward muons, divided by 2ti [|62]| : 

S(E U ) = -L|°dcos0 J d<f)exp[-z(6)/C int (Ev)}- (3-5) 



The column depth z{6) is plotted in Figure 15 of Ref. WA - In our estimates of event rates, 
we choose the interaction length that corresponds to the total (charged-current plus neutral- 
current) cross section, which is shown for v^N and v^N interactions as the solid curves in 
Figures |] and 0. The shadow factor calculated using the CTEQ3-DIS parton distributions 
is shown in Figure 20 of Ref. [|12|1 . 

The rate at which upward-going muons can be observed in a detector with effective area 
A is 

Rate = A / dE v P^E V - E™ n )S{E v )—-. (3.6) 

Je™ 1d ahi v 

The integrals in ( |3.6p are evaluated up to E max = 10 11 GeV except for AGN-SS91, for which 
the data files extend only to £ max = 10 9 - 8 GeV. 

Let us consider for illustration a detector with effective area A = 0.1 km 2 . We show in 



Tables [TTI| and [TV] the annual event rates for upward-going muons with observed energies 
exceeding 1 TeV and 10 TeV, respectively. We tabulate rates for the full upward-going solid 
angle of 2tc, as well as for the detection of "nearly horizontal" muons with nadir angle 9 
between 60° and 90°. The predicted event rates, shown here for the CTEQ4-DIS parton 
distributions, are very similar for other modern parton distributions. 

The AGN fluxes introduced in A| all yield significant rates for upward-going muon 
events, but the events induced by atmospheric neutrinos constitute an important back- 
ground, especially for E™ m = 1 TeV. The highest signal rates arise in the AGN-SS91 
model. At the lower muon energy threshold, E™ m = 1 TeV, this model yields rates about 
half the expected background from the interactions of atmospheric neutrinos. The expected 
background exceeds the signal rates for other models with the muon energy threshold set at 
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1 TeV, but fades as the threshold is raised to 10 TeV. When fully deployed, neutrino ob- 
servatories of the current generation will have instrumented areas A « 0.02 km 2 . Although 
they should register a handful of AGN neutrinos, it will be difficult to distinguish them from 
the atmospheric background. The energy dependence of the upward-going muon rate will 
be an important discriminant for separating atmospheric and extraterrestrial sources. 

The GRB-WB flux yields a small number of upward-going muon events per year in a 
0.1-km 2 detector. To pick these out from the background, it will be essential to correlate 
them in time and position with gamma-ray bursts. The TD-SLSC flux appears unobservable 
in a 0.1-km 2 detector. 

The ratio of the "nearly horizontal" rate (60° < 9 < 90°) to the full upward rate charac- 
terizes the attenuation of the incident neutrino flux in the Earth. The energy dependence of 
the neutrino interaction length shown in Figures ^| and ^ is reflected in the angular depen- 
dence of the shadow factor, for given neutrino energy: the greater the neutrino energy, the 
more the incident neutrino flux is attenuated in its passage through the Earth. We recall 
that at 9 = 60°, the column depth is about 20% of the vertical column depth [fL2| . The ATM 
and AGN-M95 fluxes fall rapidly with neutrino energy, so rates are dominated by neutrino 
energies near the muon energy threshold, and attenuation is not severe for E™ m = 1 TeV. 
For these model fluxes, the upper half of the 2tt solid angle for upward events contributes 
between 52% and 58% of the total number of upward-going events. The other model fluxes 
decrease more slowly with neutrino energy. Accordingly, the effect of shadowing on the an- 
gular distribution is more pronounced. For the AGN-SS91 and AGN-M95 fluxes, the "nearly 
horizontal" solid angle contributes 76% and 87% of the total. The relatively stiff GRB-WB 
flux also entails significant shadowing: 67% of the events come from the Itt solid angle just 
below the horizon. 

The importance of shadowing increases as the muon energy threshold is raised. For 
£™ n = 10 TeV, between 59% and 89% of the upward rate comes from the upper half of the 
upward solid angle. 

If we raise the muon energy threshold to 100 TeV, the atmospheric neutrino background 
is essentially eliminated. At the same time, the projected signal rates, shown in Table 
for a detector with A — 0.1 km 2 , reflect the influence of shadowing in the Earth in two 
ways. First, the angular distributions favor shallow angles still more strongly. Between 
77% and 94% of the extraterrestrial neutrino signal comes from the "nearly horizontal" 
wedge. Second, the shadow factor S(E U ) drops rapidly above 100 TeV. It is equal to 0.64, 
0.34, 0.16, 0.071 for E v = 100 TeV, 1 PeV, 10 PeV, and 100 PeV. Even so, the AGN-SS91 
and AGN-P96 upward rates should be observable in a detector with an effective ares of 
A = 0.1 km 2 . 

The penalty of shadowing overcomes the advantages of looking at upward-going muons as 
the muon-energy threshold is raised. For very high thresholds, it will be necessary to observe 
downward-going muons produced by interactions within the instrumented volume. We show 
in Table [VI| the annual downward event rates for v^N and v^N charged-current interactions 
in a 1-km 3 volume of water. The rates are encouraging, provided that downward-going 
events can be observed efficiently. 

We discussed the effect of different choices for the parton distribution functions in 
GQRS96 for energy thresholds of 1 and 10 TeV. We concluded that the upward muon 
event rate is insensitive to the choice of modern (post-HERA) parton distribution. We have 
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evaluated the upward muon event rates for E™ m = 100 TeV with the MRS D_' |63[] and 
CTEQ3-DLA |64|] parton distributions. The event rates are essentially equal to those shown 
in Table [V|. 

We indicated in Section [II C| that all standard parton distribution functions yield very 
similar cross sections up to E u m 10 7 GeV, so only with a muon threshold E™ m > 10 7 GeV 
might one distinguish between parton distribution functions. The angular distribution of 
the upward muons is a measure of attenuation as a function of column depth, and thus of 
the cross section. We show in Figure |10| the differential shadow factor (1/2tt) dS/dcos6 
for three neutrino energies and for three sets of parton distributions. At E v = 10 3 GeV, 
there is no appreciable shadowing at any angle. By E v = 10 6 GeV, S(E U ) 0.3, and the 
nearly vertical (upward-going) events are depleted in comparison to nearly horizontal events. 
However, the angular distribution of events will be quite similar for our nominal set of parton 
distributions (CTEQ4-DIS) and for the CTEQ3-DLA and MRS-D_' distributions. At still 
higher energies, events will be observed only near the horizon. The angular distribution 
is steepest for the MRS-D_' distributions, which yield the largest vN cross sections, and 
shallowest for the CTEQ3-DLA distributions, which yield the smallest vN cross sections. 
Given the low rates we anticipate for E v > 10 8 GeV, a discriminating measurement would 
demand a prohibitively large instrumented volume. 



C. Neutrino-electron scattering 

Observations of electron neutrino interactions at higher energies and large target volumes 
may yield insights into the ultrahigh energy neutrino flux and the high energy vN cross 
section. Generically, the v e + v e flux from the 7r — > fi — > e decay chain is one-half the 
+ flux. For the rate estimates presented here, we use the fluxes of Figure || multiplied 
by 0.5. As a cautionary note, we call attention to the contention of Rachen and Meszaros 
iH that muon cooling within astrophysical sources may reduce the UHE + flux by a 



factor of two and effectively eliminate the UHE v e + 9 e flux. In this subsection, we shall first 
evaluate rates for resonant W production. Then we consider the possibility of observing 
nearly horizontal air showers induced by neutrino interactions in the atmosphere. 

Resonant v e e — > W~ production occurs for E D ~ 6.3 PeV. In an effective volume V e s, 
the contained event rate is 

10 r {M w +2T w f/2m dN D 

Rate = - N A V eS dE^Jft,.) — ^ . (3.7) 

18 J(M w -2Y w ) 2 /2m dEp e 

Downward resonant v e e — > W~ rates are shown in Table |V1J| for an effective volume of 1 km 3 . 
To assess potential backgrounds to the detection of v e e W~ — > [iv^ and v e e — > W~ — > 
hadrons, we also show the downward (and upward) (z/ M , v^N charged- current and neutral- 
current events that occur for E v > 3 PeV. The v e N downward and upward interaction 
rates and the v e N downward interaction rates may be obtained from the {y^, rates 
by scaling the incident fluxes. The v e N upward interaction rates are reduced by the short 
v e e interaction length near resonance. All of these backgrounds to the identification of 
v e e — > W~ — > or hadrons are themselves evidence for extraterrestrial neutrinos. We 
conclude, as in GQRS96, that resonant W~ production will be difficult to extract from 



13 



the neutrino-nucleon interaction background. The short v e interaction length at resonance 
(cf. Figure ||) means that the flux of electron antineutrinos is extinguished for neutrinos 
traversing the Earth. 

Let us now take up the observability of neutrino interactions in the atmosphere. A 
neutrino normally incident on a surface detector passes through a column of density of 
1 033 cmwe, while a neutrino arriving along the horizon encounters a column of about 
36 000 cmwe. Both amounts of matter are orders of magnitude smaller than the neutrino 
interaction lengths summarized in Figures so the atmosphere is essentially transparent 
to neutrinos. 

However, the horizontal path length low in the atmosphere is not tiny compared with the 
depth available for the production of contained events in a water or ice Cerenkov detector 
|66| , so it is worth asking what capabilities a large-area air-shower array might have for the 



study of UHE neutrino interactions. The proposed Pierre Auger Cosmic Ray Observatory 
P^| , which would consist of an array of water Cerenkov tanks dispersed over a large land 
area, is designed to detect showers of particles produced in the atmosphere. Proton- and 
photon-induced showers are typically produced high in the atmosphere. Nearly horizontal 
events with shower maxima near the surface array are more likely to arise from neutrino 
interactions than from p-Air or 7- Air collisions [B7|. 



The acceptance A of the Auger ground array, which has dimensions of volume times solid 
angle, has been evaluated by several authors p4||68| , |69l . We adopt the Billoir's estimate |34 



as shown in Figure 11, to compute the event rate 



Rate = iV APair / dE sh d y-f^(E u ,y)A(E sh ) . (3.8) 
JE th Jo dE u ay 



The values of _E max are the same as for our calculation of upward muon rates in § pil B[ For 
(z/ e + z/ e )iV charged- current interactions, we take the shower energy to be the sum of hadronic 
and electromagnetic energies, E^ = E v . For (z^ + z^JiV charged- current interactions and for 
neutral-current interactions, we take the shower energy to be the hadronic energy, E^ = yE v . 
The resulting event rates, calculated using our canonical CTEQ4-DIS parton distributions, 
are shown in Table |V111] for (z/ e + v e )N neutral-current interactions and for iv^ + v^N 
charged- current interactions, for two shower thresholds. The iv^ + v^)N neutral-current 
rates are twice those shown for the (z/ e + v e )N case. In Table |X| we show our evaluation of 
the (z/ e + v e )N charged-current rates for three different sets of parton distributions. 

The (u e +u e )N neutral-current event rates are typically less than 15% of the corresponding 
charged-current rates, reflecting a combination of smaller cross sections and a falling flux. 
The (z/ M + Vy) N charged-current rates are a factor of ~ 4 larger than the (z/ e + v e )N neutral- 
current rates, because of the larger flux and larger cross section. The inelasticity parameters 
{y) are approximately equal for neutral-current and charged-current interactions. 

The largest rates for neutrino- induced horizontal air showers arise from [v e + u e )N 
charged-current interactions, for which E^ w E v . In one year, a few to tens of horizon- 
tal (z/ e + z/ e )iV —>■ e T + anything events may be observed in the Auger detectors, assuming 
the modern estimates of AGN neutrino fluxes. The AGN-SS91, GRB-WB, TD-SLSC, and 
TD-WMB16 fluxes yield fractions of an event per year, while the TD-WMB12 flux yields 
an event or two. 
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Given the high thresholds that must be set for detection, the expected event rates are 
dependent on the choice of parton distribution functions. The D_' rates are approximately 
a factor of two larger than the CTEQ3-DLA rates. If the absolute normalization and en- 
ergy behavior of the AGN fluxes could be established in underground experiments at lower 
energies, the Auger experiment might suggest distinctions among the various high-energy 
extrapolations of the cross sections. 

IV. SUMMARY AND OUTLOOK 

We have presented new calculations of the cross sections for neutrino-nucleon charged- 
current and neutral-current interactions. The new cross sections are at most 25% smaller 
than those of GQRS96, with the deviation largest at the highest energy considered here, 
10 21 eV. By varying the extrapolations of the small- x behavior of the parton distribution 
functions, we find that the uncertainty in the vN cross section is at most a factor of 2 ±1 at 
the highest energies. All modern sets of parton distribution functions give comparable cross 
sections for energies up to 10 16 eV. 

We have estimated event rates for several energy thresholds and detection methods, 
using a variety of models for the neutrino fluxes from AGNs, gamma-ray bursters, topological 
defects, and cosmic-ray interactions in the atmosphere. In v^N — > fiX interactions, requiring 
a muon energy above 10 TeV reduces the atmospheric background enough to permit the 
observation of upward-going muons for the AGN-SS91 and AGN-P96 fluxes. These models 
yield tens to hundreds of events per year for detectors of 0.1 km 2 effective area. The GRB- 
WB flux emerges at a higher threshold, but suffers from a small event rate. 

Event rates for downward muons above 100 TeV from neutrinos are substantial in 1 km 3 , 
except for the TD models. Resonant W boson production will be difficult to distinguish from 
the vN interaction background. For the Pierre Auger Cosmic Ray Observatory, the most 
promising rates arise from (u e , u e )N charged- current interactions in the AGN-M95 and 
AGN-P96 models. By combining measurements of the upward-going muon rate at lower 
energies with air-shower studies at the highest energies, it may be possible to distinguish 
among alternative high-energy extrapolations of the vN cross section. 

The origins of the highest energy cosmic rays are not well understood, but cosmic rays 
should be accompanied by very high energy neutrinos in all models. The absolute normal- 
ization and energy dependence of the fluxes vary from model to model. Neutrino telescopes 
ultimately will probe extraterrestrial accelerator sources. We expect that detectors with 
effective areas on the order of 0.1 km 2 will yield significant clues to aid in our understanding 
of physics to the 10 20 -eV energy scale. 
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TABLES 



TABLE I. Charged-current and neutral-current cross sections and their sum for vN interactions 
according to the CTEQ4-DIS distributions. 



E v [GeVl 

v L J 




[cm 2 ] 


<^NC 


[cm 2 ] 


(Ttot 


[cm 2 ] 


1.0 x 10 1 


0.7988 


X 


1 n— 37 

10 31 


0.2492 


X 


1 n — 37 

10 61 


0.1048 


X 


10 36 


2.5 x 10 1 


0.1932 


X 




0.6033 


X 


1 n— 37 

10 31 


0.2535 


X 


i r\— 36 

10 30 


6.0 x 10 1 


0.4450 


X 




0.1391 


X 


10 -36 


0.5841 


X 


lQ-36 


1.0 x 10 2 


0.7221 


X 


1 n— 3fi 

10 M 


0.2261 


X 


-i r>— 3fi 

10 30 


0.9482 


X 


i rv— 3fi 

10 M 


2.5 x 10 2 


0.1728 


X 


lQ-35 


0.5430 


X 


lQ-36 


0.2271 


X 


10~ 35 


6.0 x 10 2 


0.3964 


X 




0.1255 


X 


IO -35 


0.5219 


X 


1Q-35 


1.0 x 10 3 


0.6399 


X 


i r> — 3^ 

10 M 


0.2039 


X 


i n — 3^ 

10 30 


0.8438 


X 


i r> — 3^ 

10 30 


2.5 x 10 3 


0.1472 


X 


i n — 34 

10 M 


0.4781 


X 


i n— 3^ 

10 30 


0.1950 


X 


i n — 34 

10 134 


6.0 x 10 3 


0.3096 


X 


i r,~ 34 

10 d4 


0.1035 


X 


i r> — 34 

10 d4 


0.4131 


X 


i n — 34 

10 d4 


1.0 x 10 4 


0.4617 


X 


-i n— 34 

10 * 


0.1575 


X 


i n— 34 

10 3 


0.6192 


X 


i n — 34 

10 3 


2.5 x 10 4 


0.8824 


X 


-i n— 34 

10 * 


0.3139 


X 


-i n— 34 

10 3 


0.1196 


X 


1Q-33 


6.0 x 10 4 


0.1514 


X 


-I n— 33 

10 M 


0.5615 


X 


i n — 34 

10 134 


0.2076 


X 


1 n — 33 

10 M 


1.0 x 10 5 


0.2022 


X 


i n— 33 
10 M 


0.7667 


X 


i n — 34 

10 d4 


0.2789 


X 


1 n — 33 

10 


2.5 x 10 5 


0.3255 


X 


1 n — 33 

10 M 


0.1280 


X 


i n — 33 
10 M 


0.4535 


X 


1 n— 33 

10 


6.0 x 10 5 


0.4985 


X 


10 33 


0.2017 


X 


10 -33 


0.7002 


X 


i n — 33 

10 


1.0 x 10 6 


0.6342 


X 


IO -33 


0.2600 


X 


10 -33 


0.8942 


X 


10~ 33 


2.5 x 10 6 


0.9601 


X 


i n — 33 

10 33 


0.4018 


X 


i n — 33 
10 ^ 


0.1362 


X 


i n — 32 

10 


6.0 x 10 6 


0.1412 


X 


i n — 32 

10 


0.6001 


X 


i n — 33 

10 


0.2012 


X 


i n — 32 
10 


1.0 x 10 7 


0.1749 


X 


i n — 32 

10 


0.7482 


X 


i r\— 33 
10 ^ 


0.2497 


X 


i n — 32 

10 


2.5 x 10 7 


0.2554 


X 


i n— 32 

10 3 


0.1104 


X 


i r>— 32 
10 


0.3658 


X 


i n — 32 
10 


6.0 x 10 7 


0.3630 


X 


i n — 32 

10 3 


0.1581 


X 


i n— 32 

10 3 


0.5211 


X 


i n — 32 
10 


1.0 x 10 s 


0.4436 


X 


-I n— 32 

10 6Z 


0.1939 


X 


i n — 32 
10 ^ 


0.6375 


X 


i n — 32 
10 ^ 


2.5 x 10 8 


0.6283 


X 


i r> — 32 

10 3Z 


0.2763 


X 


i n — 32 

10 


0.9046 


X 


i n — 32 

10 


6.0 x 10 8 


0.8699 


x 


1Q -32 


0.3837 


x 


1Q -32 


0.1254 


x 


10~ 31 


1.0 x 10 9 


0.1049 


X 


io- 31 


0.4641 


X 


1Q -32 


0.1513 


X 


io- 31 


2.5 x 10 9 


0.1466 


X 


10~ 31 


0.6490 


X 


1Q -32 


0.2115 


X 


10~ 31 


6.0 x 10 9 


0.2010 


X 


io- 31 


0.8931 


X 


1Q -32 


0.2903 


X 


io- 31 


1.0 x 10 10 


0.2379 


X 


io- 31 


0.1066 


X 


io- 31 


0.3445 


X 


io- 31 


2.5 x 10 10 


0.3289 


X 


io- 31 


0.1465 


X 


io- 31 


0.4754 


X 


io- 31 


6.0 x 10 10 


0.4427 


X 


io- 31 


0.1995 


X 


io- 31 


0.6422 


X 


io- 31 


1.0 x 10 11 


0.5357 


X 


io- 31 


0.2377 


X 


io- 31 


0.7734 


X 


io- 31 


2.5 x 10 11 


0.7320 


X 


io- 31 


0.3247 


X 


io- 31 


0.1057 


X 


10 -30 


6.0 x 10 11 


0.9927 


X 


10~ 31 


0.4377 


X 


io- 31 


0.1430 


X 


10 -30 


1.0 x 10 12 


0.1179 


X 


10 -30 


0.5196 


X 


10~ 31 


0.1699 


X 


10 -30 
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TABLE II. Charged-current and neutral-current cross sections and their sum for i/N interac- 
tions according to the CTEQ4-DIS distributions. 



E v [GeV] 




[cm 2 ] 




[cm 2 ] 


0"tot 


[cm 2 ] 


1.0 x 10 1 


0.3936 


X 


io- 37 


0.1381 


X 


io- 37 


0.5317 


X 


io- 37 


2.5 x 10 1 


0.9726 


X 


io- 37 


0.3403 


X 


10~ 37 


0.1313 


X 


1Q -36 


6.0 x 10 1 


0.2287 


X 


1Q -36 


0.7982 


X 


10~ 37 


0.3085 


X 


1Q -36 


1.0 x 10 2 


0.3747 


X 


1Q -36 


0.1307 


X 


1Q -36 


0.5054 


X 


10~ 36 


2.5 x 10 2 


0.9154 


X 


1Q -36 


0.3193 


X 


1Q -36 


0.1235 


X 


1Q -35 


6.0 x 10 2 


0.2153 


X 


1Q -35 


0.7531 


X 


1Q -36 


0.2906 


X 


1Q -35 


1.0 x 10 3 


0.3542 


X 


1Q -35 


0.1243 


X 


1Q -35 


0.4785 


X 


1Q -35 


2.5 x 10 3 


0.8548 


X 


1Q -35 


0.3026 


X 


1Q -35 


0.1157 


X 


1Q -34 


6.0 x 10 3 


0.1922 


X 


1Q -34 


0.6896 


X 


10 -35 


0.2612 


X 


10~ 34 


1.0 x 10 4 


0.3008 


X 


1Q -34 


0.1091 


X 


1Q -34 


0.4099 


X 


10 -34 


2.5 x 10 4 


0.6355 


X 


1Q -34 


0.2358 


X 


1Q -34 


0.8713 


X 


10~ 34 


6.0 x 10 4 


0.1199 


X 


10~ 33 


0.4570 


X 


lfJ -34 


0.1656 


X 


10~ 33 


1.0 x 10 5 


0.1683 


X 


1Q -33 


0.6515 


X 


1Q -34 


0.2334 


X 


1Q -33 


2.5 x 10 5 


0.2909 


X 


10 -33 


0.1158 


X 


1Q -33 


0.4067 


X 


1Q -33 


6.0 x 10 5 


0.4667 


X 


1Q -33 


0.1901 


X 


1Q -33 


0.6568 


X 


1Q -33 


1.0 x 10 6 


0.6051 


X 


1Q -33 


0.2493 


X 


1Q -33 


0.8544 


X 


1Q -33 


2.5 x 10 6 


0.9365 


X 


1Q -33 


0.3929 


X 


1Q -33 


0.1329 


X 


1Q -32 


6.0 x 10 6 


0.1393 


X 


1Q -32 


0.5930 


X 


1Q -33 


0.1986 


X 


1Q -32 


1.0 x 10 7 


0.1734 


X 


1Q -32 


0.7423 


X 


1Q -33 


0.2476 


X 


1Q -32 


2.5 x 10 7 


0.2542 


X 


1Q -32 


0.1100 


X 


IO" 32 


0.3642 


X 


1Q -32 


6.0 x 10 7 


0.3622 


X 


1Q -32 


0.1578 


X 


IO" 32 


0.5200 


X 


1Q -32 


1.0 x 10 s 


0.4430 


X 


1Q -32 


0.1937 


X 


1Q -32 


0.6367 


X 


1Q -32 


2.5 x 10 8 


0.6278 


X 


1Q -32 


0.2762 


X 


1Q -32 


0.9040 


X 


1Q -32 


6.0 x 10 s 


0.8696 


X 


10 -32 


0.3836 


X 


1Q -32 


0.1253 


X 


io- 31 


1.0 x 10 9 


0.1050 


X 


io- 31 


0.4641 


X 


1Q -32 


0.1514 


X 


10~ 31 


2.5 x 10 9 


0.1464 


X 


10~ 31 


0.6489 


X 


1Q -32 


0.2113 


X 


10~ 31 


fi D v If) 9 

U.U A 1U 




X 


io- 31 


u.oyoi 


X 


1Q -32 




X 


io- 31 


1.0 x 10 10 


0.2406 


X 


io- 31 


0.1066 


X 


io- 31 


0.3472 


X 


io- 31 


2.5 x 10 10 


0.3286 


X 


io- 31 


0.1465 


X 


io- 31 


0.4751 


X 


io- 31 


6.0 x 10 10 


0.4481 


X 


IO" 31 


0.1995 


X 


io~ 31 


0.6476 


X 


10~ 31 


1.0 x 10 11 


0.5335 


X 


10~ 31 


0.2377 


X 


io~ 31 


0.7712 


X 


10~ 31 


2.5 x 10 11 


0.7306 


X 


IO" 31 


0.3247 


X 


io~ 31 


0.1055 


X 


10 -30 


6.0 x 10 11 


0.9854 


X 


io~ 31 


0.4377 


X 


10~ 31 


0.1423 


X 


10 -30 


1.0 x 10 12 


0.1165 


X 


lfJ -30 


0.5195 


X 


IO" 31 


0.1685 


X 


10 -30 
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TABLE III. Upward fi + + /x - event rates per year arising from v^N and U^N interactions in 
rock, for a detector with effective area A = 0.1 km 2 and muon energy threshold E™ m = 1 TeV. The 
rates are shown integrated over all angles below the horizon and restricted to "nearly horizontal" 
nadir angles 60° < 9 < 90°. 



Flux 



nadir angular acceptance 
0° < 9 < 90° 60° < 9 < 90 c 



ATM H| 

ATM @ + charm 
AGN-SS91 || 
AGN-M95 Oy) @|] 
AGN-P96 (p7) 
GRB-WB M 



TD-SLSC H 



TD-WMB12 |5| 



TD-WMB16 H 



1100 


570 


1100 


570 


500 


380 


31 


18 


45 


39 


12 


8.1 


0.005 


0.0046 


0.50 


0.39 


0.00050 


0.00039 



TABLE IV. Upward /i + + /i - event rates per year arising from v„N and v^N interactions in 
rock, for a detector with effective area A = 0.1 km 2 and muon energy threshold E™ m = 10 TeV. 
The rates are shown integrated over all angles below the horizon and restricted to "nearly horizon- 
tal" nadir angles 60° < 9 < 90°. 









nadir anj 


^ular acceptance 


Flux 






0° < 9 < 90° 


60° < 9 < 90° 


ATM [ 


58 


1 


17 


10 


ATM [ 


58 


+ charm |pQ] 


19 


11 


AGN-SS91 || 


270 


210 


AGN-M95 Oy) @ 


5.7 


4.3 


AGN-P96 (p7) || 


28 


25 


GRB-WB @ 


5.4 


4.0 
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TABLE V. Upward event rates per year arising from v^N and VnN interactions in rock, 

for a detector with effective area A = 0.1 km 2 and muon energy threshold E™ m = 100 TeV. The 
rates are shown integrated over all angles below the horizon and restricted to "nearly horizontal" 
nadir angles 60° < 6 < 90°. 



Flux 



0° < 9 < 90° 



nadir angular acceptance 



60° < 6 < 90 c 



ATM H| 

ATM @ + charm 
AGN-SS91 || 
AGN-M95 Oy) @|] 
AGN-P96 (p7) 
GRB-WB M 



0.13 

0.21 
85 

1.6 
13 

1.2 



0.09 

0.16 
73 

1.5 
12 

1.0 



TABLE VI. Downward ^ 
1 km 3 of water. 



// events per year arising from v^N and v^N interactions in 







Muon-energy threshold, E™ m 




Flux 


100 TeV 


1 PeV 


3 PeV 



ATM p| 



ATM |58| + charm 
AGN-SS91 || 
AGN-M95 (p7) || 
AGN-P96 (p7) " 
GRB-WB @ 
TD-SLSC p[ 



TD-WMB12 H 



TD-WMB16 g 



0.85 

2.6 
520 
16 
100 

7.7 

0.037 

1.1 

0.00087 



0.0054 
0.050 
120 

11 

50 
1.9 
0.032 
0.74 
0.00050 



0.00047 

0.0076 
42 

8.7 
31 

0.93 

0.029 

0.58 

0.00035 
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TABLE VII. Downward resonant v e e — ► W~ events per steradian per year for a detector with 
effective volume V e R = 1 km 3 . Also shown are the downward (upward) potential background rates 
from v^N and VuN interactions induced by neutrinos with E v > 3 PeV. 



Flux 


z/ e e -> z/^// 


f e e — ► hadrons 




)N CC 




)N NC 


AGN-SS91 || 


6 


41 


29 


(5.2) 


13 


(2.3) 


AGN-M95 Oy) |j 


| 0.1 


0.6 


2.3 


(0.21) 


1.1 


(0.095) 


AGN-P96 (p7) [48 


| 1.2 


7.8 


12 


(1.6) 


5.2 


(0.69) 


GRB-WB [0 


0.06 


0.4 


0.43 


(0.065) 


0.19 


(0.029) 


TD-SLSC |51] 


0.001 


0.0074 


0.0059 


(0.00031) 


0.0028 


(0.00014) 



TABLE VIII. Annual neutral-current {y e , v e )N and charged-current (u„, Vn 
the Pierre Auger Cosmic Ray Observatory. 



)N event rates for 



Flux 


E sh > 10 8 GeV 


E sh > 10 9 GeV 




(u e , u e )N NC (zv v^N CC 


(u e , u e )N NC z^)iV CC 



AGN-SS91 
AGN-M95 (p7) 
AGN-P96 (p7) 
GRB-WB lEOI 



TD-SLSC H 
TD-WMB12 H 
TD-WMB16 pi 



0.0045 

0.65 

0.74 

0.038 

0.013 

0.15 

0.000026 



0.019 

2.7 

3.1 

0.16 

0.052 

0.59 

0.00011 



0.000006 

0.26 

0.13 

0.020 

0.010 

0.11 

0.000011 



0.000024 

1.1 

0.53 

0.085 

0.042 

0.44 

0.000046 
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TABLE IX. Annual event rates in the Pierre Auger Cosmic Ray Observatory for horizontal air 
showers induced by (v e ,D e )N charged-current interactions. 



Flux 




Parton distributions 


E sh > 10 8 GeV 


E sh > 10 9 GeV 


AGN-SS91 |j 


I 


CTEQ4-DIS 

G1EQ3-DLA 

D_' 


0.15 
0.13 
0.23 


0.00026 

r\ r\f\r\c\c\ 

0.00022 
0.00051 


AGN-M95 (jry) @ 


CTEQ4-DIS 

/~ir"nT7i/"\ O TM A 

G1EQ3-DLA 
D_' 


6.1 
5.3 
12 


3.3 
2.8 
7.5 


AGN-P96 (p7) H 


CTEQ4-DIS 

CTEQ3-DLA 

D_' 


8.9 
7.9 

16 


2.6 
2.2 
5.4 


GRB-WB |H 




CTEQ4-DIS 

CTEQ3-DLA 

D_' 


0.31 
0.27 
0.67 


0.18 
0.16 
0.45 


TD-SLSC |51 




CTEQ4-DIS 

CTEQ3-DLA 

D_' 


0.068 
0.056 
0.18 


0.061 
0.051 
0.17 


TD-WMB12 


52] 


CTEQ4-DIS 

CTEQ3-DLA 

D_' 


0.85 
0.72 
2.1 


0.71 
0.60 
1.9 


TD-WMB16 


52] 


CTEQ4-DIS 

CTEQ3-DLA 

D_' 


0.00024 
0.00021 
0.00049 


0.00014 
0.00012 
0.00032 
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FIG. 1. Cross sections for vgN interactions at high energies, according to the CTEQ4-DIS 
parton distributions: dashed line, a(viN — > U£ + anything); thin line, a{y^N — ► l~ + anything); 
thick line, total (charged-current plus neutral-current) cross section. 
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FIG. 2. Ratio of the charged-current cross section shown in Figure |l|, calculated using the 
CTEQ4-DIS parton distributions, to the charged-current cross section of Ref. |l2| calculated using 
the CTEQ3-DIS parton distributions. 
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E„ [GeV] 

FIG. 3. Cross sections for v^N interactions at high energies, according to the CTEQ4-DIS 
parton distributions: dashed line, a{v^N — ► + anything); thin line, a(vgN — ► l + + anything); 
thick line, total (charged-current plus neutral-current) cross section. 
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FIG. 4. Cross sections for v^N interactions at high energies, according to the CTEQ4-HJ 
parton distributions: dashed line, a{v^N — > v g + anything); thin line, a{yiN — ► l~ + anything); 
thick line, total (charged-current plus neutral-current) cross section. 
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FIG. 5. Ratio of the charged-current cross section shown in Figure calculated using the 
CTEQ4-HJ parton distributions, to the charged-current cross section of Figure [j] calculated using 
the CTEQ4-DIS parton distributions. 
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E„ [GeV] 

FIG. 6. Interaction lengths for neutrino interactions on nucleon targets: dotted line, 
charged-current interaction length; dashed line, neutral-current interaction length; solid line, total 
interaction length, all computed with the CTEQ4-DIS parton distributions. The dot-dashed curve 
shows the charged-current interaction length based on the EHLQ structure functions with Q 2 held 
fixed at Ql = 5 GeV 2 . 




10 100 1000 10 4 10 5 10 6 10 7 10 8 10 9 10 10 10 n 10 



E v [GeV] 

FIG. 7. Interaction lengths for antineutrino interactions on nucleon targets: dotted line, 
charged-current interaction length; dashed line, neutral-current interaction length; solid line, total 
interaction length, all computed with the CTEQ4-DIS parton distributions. 
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E y [GeV] 

FIG. 8. Interaction lengths for neutrino interactions on electron targets. At low energies, from 
smallest to largest interaction length, the processes are (i) v e e — > hadrons, (ii) i/„e — >• iiv e , (hi) 
z^ e e -> z^ e e, (iv) v e e -> z7„/i, (v) ^ e e -> P e e, (vi) z/„e -> f„e, (vii) P^e -> z?„e. [From Ref. @] 




i:, [GeV] 

FIG. 9. Muon neutrino plus antineutrino fluxes scaled by neutrino energy at the Earth's surface. 
Solid lines represent AGN fluxes. In decreasing magnitude at E u = 10 3 GeV, they are AGN-M95, 
AGN-SS91 scaled by 0.3, and AGN-P96 (jry). The dashed lines, in the same order, represent the 
GRB-WB, TD-WMB12, TD-WMB16, and TD-SLSC fluxes. The dotted line is the angle-averaged 
atmospheric (ATM) neutrino flux. 



31 




FIG. 10. Differential shadow factor (l/2n) dS (E u ) / d cos 9 versus nadir angle 9, for 
E v = 10 3 GeV (solid line), 10 6 GeV (dot-dashed lines), and 10 8 GeV (dashed lines). For each 
neutrino energy, the thick line correspondes to the CTEQ4-DIS parton distributions; the upper 
and lower satellite lines correspond to the CTEQ3-DLA and MRS-D_' parton distributions, re- 
spectively. 




FIG. 11. Acceptance versus shower energy for nearly horizontal air showers in the proposed 
Pierre Auger Cosmic Ray Observatory. The solid line for E s ^ < 10 10 GeV shows the acceptance 



evaluated by Billoir [34]; the dashed line is our projection to higher energies. 
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